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Abstract 

This paper presents a formulation of Noether's theorem for frac- 
tional classical fields. We extend the variational formulations for frac- 
tional discrete systems to fractional field systems. By applying the 
variational principle to a fractional action S, we obtain the fractional 
Euler-Lagrange equations of motion. Considerations of the Noether's 
variational problem for discrete systems whose action is invariant un- 
der gauge transformations will be extended to fractional variational 
problems for classical fields. The conservation laws associated with 
fractional classical fields are derived. As an example we present the 
conservation laws for the fractional Dirac fields. 

Keywords: Fractional calculus; fractional variational principle; Noether's 
theorem. 
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1 Introduction 

Fractional calculus generalized the classical calculus and it has many impor- 
tant applications in various fields of science and engineering [TJ El El SI El E] . 



1 On leave of absence from Al-Azhar University-Gaza, Email: smuslih@ictp.it 
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Also, it plays an important rule in the understanding of both conservative 
and non conservative systems [7J. Particularlythe fractional and anomalous 
dynamics has experienced to be a framework in the theory of stochastic pro- 
cesses j8] . The physical and geometrical meaning of the fractional derivatives 
has been investigated by several authors [HI [10]. It is know what is the physical 
interpretation of the Stieltjes integral, namely the Stieltjes integral can be in- 
terpreted as the real distance passed by a moving object. Within this physical 
interpretation the Riemann-Liouville and Caputo have the same interpreta- 
tion and therefore can be used successfully , for example, within fractional 
variational principles [H1II21II2H0H3II3II3IIE]. On the other hand, by 
making use of examples from viscoelasticity it was shown that it is possible 
to attribute physical meaning to initial conditions within Riemann-Liouville 
fractional derivatives [19]. The fractional operators are particular cases of 
non-local operators and during the last decade several interesting research 
was done in this direction [20] |2T| |22| 123] . 

Conservation laws play an important role in theoretical physics. The 
conservation of energy, momentum, and angular momentum are fundamen- 
tal laws that any theory has to guarantee if it is to give a valid description 
of nature. Besides, physical systems often possess further conserved quan- 
tities, such as charge, isospin. Noether's theorem [23] is a central result of 
the calculus of variation and explain all the conservation laws based upon 
the action principle. It is a very general result, asserting that "the conserva- 
tion laws are natural consequences of the symmetry properties of a system. 
For each continues transformation of the coordinates and/or the fields under 
which "the action is invariant" the existence of a conserved quantity can be 
deduced". For example, the conservation of energy, momentum, and angu- 
lar momentum are based on the invariance of the action S under temporal 
and spatial translations and under rotation in space. Similarly the conser- 
vation of charge follows from an invariance of the action S under the phase 
transformations. 

Following the above progress, one can naturally ask: " Can we apply the 
Noether's variational problem to Lagrangian systems with fractional deriva- 
tives." This is exactly what is done in [23] |2"5] . Frederic and Torres 12"6] 
examine the Noether's variational problem to Lagrangian systems with finite 
degrees of freedom. At this stage, one may ask, " Can we extend the Noether's 
variational problem to fractional classical fields". To address this, we first 
define a variational problem for classical fields [18]. A variational problem is 
a problem that requires finding the extremum of a functional which may be 
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subjected to algebraic and/or dynamic constraints. If either the functional 
and/or the algebraic/dynamic constraint (s) contain at least one fractional 
derivative term, then the problem is called a fractional variational problem. 
The Fractional Variational Calculus is an extension of the ordinary Varia- 
tional Calculus that deals with finding the solution of a fractional variational 
problem HH H21 H31 [HI H51 HQ HZ]- We will define Noether's variational 
problem to fractional Lagrangian systems with infinite degrees of freedom. 

An important point to be specified here is that, several definitions have 
been proposed of a fractional derivative, among those Riemann-Liouville and 
Caputo fractional derivatives are most popular. The differential equations 
defined in terms of Riemann-Liouville derivatives require fractional initial 
conditions whereas the differential equations defined in terms of Caputo 
derivatives require regular boundary conditions. For this reason, Caputo 
fractional derivatives are popular among scientists and engineers. Accord- 
ingly, we shall also develop our formulations in terms of Caputo fractional 
derivatives. However, most of the approach will also be applicable to prob- 
lems defined using other derivatives. 

The plan of this paper is as follows: In section 2, we present the Euler- 
Lagrange equations for a fractional field. Section 3 presents the Noether's 
variational theorem for fractional classical fields. In section 4, we define 
a fractional Lagrangian density function and use the theories developed in 
sections 2 and 3 to derive the conservation laws for fractional Dirac fields. 
Finally, section 5 presents conclusions. 

2 Lagrangian formulation of field systems with 
Caputo fractional derivatives 

Consider a function / depending on n variables, X\ , ■ ■ ■ , x n defined over the 
domain Q = [ai, b\] x ■ ■ ■ x [a n , b n }. Following the convention used in physics, 
we defined the left and the right partial Riemann-Liouville and Caputo frac- 
tional derivatives of order a k} < a k < 1 with respect to Xfc cLS I lol 



(+37) (*) 



r(i-a fc ) 



1 



dx k 



■X 



k fjxu ■ ■ ■ ,X k -i,U,X k+ i, ■ ■ ■ 
(x k - u) ak 




du, 



(1) 



(-37) (*) 



r(i-a fc ) 



-1 



dx k 



' 6fc fjxi, ■ ■ ■ ,x k -i,u,x k+1 , ■ ■ ■ 
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(Caat\/ m \_ 1 f Xk d u f(Xi, ■ ■ ■ , Xk-l, U, • • • , X n ) 



r(l - a k ) Ja k (x k - u)l 

and 

fCo>a.£\t \ ~ 1 f bk dufyXi, ■ ■ ■ , Xk—i, u, Xk+i, • ■ ■ , x n ) 

{-Ouj){x) = — r / t - 0.U, (4) 

V k ' V ' r(l-a fc ) Jx k (u-x k )°"> V ; 

where dx k g is the partial derivatives of g with respect to the variable x k - 
Here, in +d%, -d%, and c _d%, the meaning of various subscripts and 

superscripts need to be made clear. The subscript k and the superscript 
a indicate that the derivative is taken with respect to the variable x k and 
it is of order a k (note that we write only a for a k , and the subscript k to 
d also represents the subscript to a), the subscripts + and — prior to the 
symbol d represent the left and the right fractional derivatives respectively, 
and accordingly the limits of integrations are taken as [a k ,x k ] and [x k ,b k ]. 
Further, no superscript and the superscript C prior to the symbol d repre- 
sent the Riemann-Liouville fractional derivative and the Caputo fractional 
derivative, respectively. Superscript a is necessary here as a reminder that 
the operator d a represents a fractional derivative. When a is equal to 1, 
the superscript a can be neglected. Although, our aim in this section is to 
present the action principle for systems defined in terms of Caputo fractional 
derivatives, the Riemann-Liouville fractional derivatives are also defined here 
because they naturally arise in the formulation. 

To develop the action principle for field systems described in terms of 
fractional derivatives, consider a Lagrangian density C depends on iV distinct 
fields <j) r , (r = 1,...,N), and their fractional derivatives written as C = 
£ (4>r(%k), (+<9fc )(j) r (x k ), ( c ld k 3 )(j) r (xk),Xk)- The functional S((p) is defined as 

5(0)= / C(M^h),(+9^)Mxk)A-^)Mxk),x k )(dx k ). (5) 

Here we have used Goldstein's [27] notation. Accordingly, x k represents n 
variables x\ to x n , <p r (x k ) = (j) r (xi, ■■■ ,x n ), £(*, %d%, *, *) = £(*, ^df, ■ ■■ , 
^9", *, *), (dx k ) = dx\ ■ ■ ■ dx n , and the integration is taken over the entire 
domain Q. Other terms are defined accordingly. 

To find the necessary condition for extremum of the action functional 
defined above, consider a one parameter family of possible functions 4> r (x k ; e) 
as follows, 

4> r (x k ; e) = 4> r {x k \ 0) + erj r (x k ), (6) 



4 



where (p r (x k ;0) are the correct functions which satisfy the Hamilton's prin- 
ciple for the fractional system, r] r (x k ) are well-behaved functions that vanish 
at the endpoints, and e is an arbitrary parameter. Note that S[(p r (x k ; e)] is 
extremum at e = 0. Substituting Eq. (JSJ) into Eq. (JSJ), differentiating the 
resulting expression with respect to e, and then setting e to 0, we obtain, 

(7) 

Finally, using the formula for integration by part [12], the fact that r](x k ) is 
zero at the boundary, and a lemma from Calculus of Variations, we obtain 

Equation (J5]) is the Euler-Lagrange equation for the fractional field system. 
For a k ,(3 k — > 1, Eq. (J5]) gives the usual Euler-Lagrange equations for classical 
fields. 

It is worthwhile to mention that the same approach can be used to find 
the Euler-Lagrange equations for functionals defined in terms of Riemann- 
Liouville or mixed (Caputo and Riemann-Liouville) derivatives [28J. 



3 Noether's variational problem for fractional 
classical fields 

Like discrete systems and integer order fields, it is possible to develop Noether's 
variational problem for fields defined in terms of fractional derivatives. Let 
us consider the following space-time and fields transformations respectively 

as 

x' k = x k + 5x k , (9) 

(f>' r (x' k ) = (f)r{x k ) + S<p r (x k ). (10) 

Under the transformation (9) and (10), the resulting change in the La- 
grangian density C (c/) r (x k ), (^)^ r (x fc ), (^d^.)(/) r (x k ),x k J will be 

a 0;(4), (?W r K), (%Vr(4),4) 
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Now let us study the consequences that follow if the transformations (9) 
and (10) leave the action integral invariant, i.e., we demand 

6s(4>) = f Qf £ (^(4),(S^)0;(4),(^)0;(4),4) W) - 

J n JC((f) r (x k ), (?^)0 r (x fc ), { C _d p k )(t> r {x k ),x k ) (dx k ) = 0, (11) 

here fl' denotes the volume of integration expressed in terms of the new 
coordinates x' k . This will introduce a Jacobi determinant in the first order 

as 

d(x' k ) 



dx' k = 



d(x k ) ^ (12) 

The conservations laws in (10) are related to the symmetry transformations 
(9,10). However, in the case of gauge transformations in field theory , where 
the field have internal symmetry, we are concerned with transformations of 
the fields only (i.e., the dependent variables (10)), and not transformations 
of the space-time coordinates (the independent variables (9)) and hence we 
ignore the terms in 8x k . In this case, we have (j)' r (xk) = 4> r (%k) + &<t>r{ x k) and 
(11) becomes 



N fir fir fir 



Using the Euler Lagrange equations of motion 

in (13), we obtain the conserved quantity for fractional classical field under 
the internal symmetry transformations as 

Defining the new operators and as left and right operators 

respectively as 

C + Dt% g] = f(J^g) - (_d£f)g, (16) 

c Df[f, g] = f g ) - ( + af f) g , (i7) 
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then the conserved quantity (15) can be written as 



E C + D t 



k,r 



dC 



,$4>r 



dC 



5<fi r 



= 0. 



(18) 



For a Lagrangian density C = C (<f) r (xk), (+d k *)4> r (xk), x^j and in the limit 
a — > 1, we obtain 

DC 



E^ 



k.r 



d{d k ct)r) 



,6<f>r 



0. 



(19) 



Now, consider the following infinitesimal internal symmetry transforma- 
tions: 



<t>' r {x k ) = 4> r {x k ) + ieJ2^rs<Ps(Xk), 

s 

5x k = 0, 

the Noether's conserved quantity (15) reads as 

dC \ 



(20) 
(21) 
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ie E ■ 

k,r,s 



dC 



d{%dt<j>)) 

dC \ 



A, 



A, 



+ 



0. (22) 



As an application of this concept, let us study a complex field described by 
a Lagrangian density which is invariant under the phase transformations 



0' = + i€(f), 
(f)'* = 0* - 

This leads to the conserved quantity 

dC 



(23) 
(24) 



ie E 

k 

k 

ie E 

k 

« e E 

k 



dC 



-0\ 



dC 



k 5(? 
dC 



d{ c + d%<t>*) 
dC 



( c ld? 



em - 



-0: 



dC 



dC \ 



+ 



fc a(?af0*), 



= o. 



(25) 
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4 Noether's theorem for Dirac's field with 
Caputo fractional derivatives 

We will now apply Noether's method developed in sections 2 and 3 to derive 
fractional conservation laws for fractional Dirac field [29]. We shall limit our 
discussion to four dimensional system (the first three for space, Xi, x 2 and 
£3, and the fourth for time, £4 = it, where % is the imaginary unit. Note 
that we are considering the units that takes the speed of light equal to 1). 
The Greek indices //, A, v etc. will range from 1 to 4, the Roman indices i, 
j, k etc. will range from 1 to 3, and unless specifically stated, the repeated 
indices will represent summation. Following this convention, we propose the 
following Lagrangian density field [18] 

C = m a M + (J8£tf) (26) 

where \I> is a wave function (the bold character \l/ suggests that it could be a 
vector function. In literature, it is also called a spinor), m is the mass of the 
particle, a (= cti = a 2 = 03 = 0.4) is the order of the fractional derivative, 7 M 
are matrices, ^ is the adjoint wave function, The dimensions of ^ and 7^ 
depend on the order of the derivatives and the way the algebra of matrices 
7 M is developed [T8] . 

Using Eqs. (|HJ) and (1261) . the Euler-Lagrange equations for variables ^ 
and \I/ are given as 

7 M (+d^) + m a ^ = 0, (27) 

and 

7 M + m Q ^ = 0. (28) 

The Lagrangian density (26) is invariant under infinitesimal internal symme- 
try transformations: 

$ = ^ e ie , (29) 

^ _». = f e _ie , (30) 

where is a constant. Using Eq.(25), the conserved quantity for fractional 
Dirac's field is given by 

^rf - _d£(^7 M )^ = 0. (31) 
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(32) 



then the Dirac's conserved quantity can be written as 



(33) 



It is worthwhile to mentio that for a — > 1, then 



D% g] = 



(34) 



In this case the Dirac's continuity equation reads as 




(35) 



where 



(36) 



is the familiar electric 4-current. 

5 Conclusions 

In this paper we have presented the formulation of Noether's theorem for frac- 
tional classical fields. We used the fractional variational principle to derive 
the equations of motion for fractional classical field as well as the conserva- 
tion laws associated with the internal symmetry of gauge transformations of 
classical fields. 

As an example we have derived the continuity equation for fractional 
Dirac's field of order a. Under the limit that a — > 1, we obtain the conserved 
4-current for the regular Dirac's field. 
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